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We derive target mass corrections (TMC) for the spin-dependent nucleon structure function g\ 
and polarization asymmetry A\ in collinear factorization at leading twist. The TMCs are found 
to be significant for g\ at large xb, even at relatively high Q 2 values, but largely cancel in ^4i. A 
comparison of TMCs obtained from collinear factorization and from the operator product expansion 
shows that at low Q 2 the corrections drive the proton Ai in opposite directions. 



I. INTRODUCTION 

Understanding the transition from the perturbative 
to the nonperturbative regimes of Quantum Chromody- 
namics (QCD) remains one of the most challenging prob- 
lems in nuclear and hadron physics. Recent progress in 
describing this transition has focused on quark-hadron 
duality, which relates observables computed from quark 
and gluon degrees of freedom to those parametrized in 
terms of hadronic variables. A classic example of this is 
the phenomenon of Bloom-Gilman duality [lj, in which 
the inclusive structure functions in the region dominated 
by low-lying nucleon resonances follow deep inelastic 
structure functions describing high energy data, to which 
the resonance structure functions average Q. 

In QCD the observation of this duality can be formu- 
lated within the operator product expansion (OPE), in 
which moments of structure functions are expanded in 
inverse powers of Q 2 , the four-momentum squared of the 
exchanged photon. The leading term is given by ma- 
trix elements of twist-two local operators, and is asso- 
ciated with single parton scattering, while the 0{1/Q 2 ) 
and higher terms are related to higher twist nonpertur- 
bative multi-parton correlations. (The twist of a local 
operator in the OPE is defined as its mass dimension 
minus its spin.) The magnitude of the higher twist con- 
tributions then determines the degree to which duality 
holds 0. 

In order to reliably extract information on the duality- 
violating higher twist contributions to structure func- 
tions, it is vital to remove from the data kinematical 
corrections associated with nonzero values of Q 2 jv 2 — 
Ax 2 B M 2 /Q 2 , where v is the energy transfer, M is the 
nucleon mass, and xb = Q 2 /2Mv is the Bjorken scal- 
ing variable. While formally related to twist-two opera- 
tors [1] , these "target mass corrections" (TMCs) are sup- 
pressed by powers of M 2 /Q 2 , hence TMCs are sometimes 
inaccurately referred to as "kinematical higher twists" , 
and can obscure information on genuine higher twist 
terms. 

The importance of TMCs has been highlighted re- 
cently by high-precision data from Jefferson Lab on both 
spin-averaged and spin-dependent structure functions Q 



taken at moderate Q 2 values, Q 2 ~ 1 — 5 GeV 2 , and 
at large xb, where TMCs are most significant. Further- 
more, with high-intensity neutrino-nucleus scattering ex- 
periments planned in similar kinematics [6[, TMCs for 
weak interactions also need to be understood. 

Target mass corrections for spin-averaged nucleon 
structure functions were first considered by Georgi & 
Politzer within the OPE 0, 0], and later extended to 
the full set of electroweak structure functions [1, Q . For 
spin-dependent scattering, these were evaluated within 
the same OPE formalism in Refs. (IcI JTll|. extended to 
the full set by Blumlein & Tkabladze [g], and computed 
by Detmold [12| for the deuteron. 

One of the limitations of the OPE formulation of TMCs 
is the so-called "threshold problem", in which the tar- 
get mass corrected structure functions remain nonzero 
at xb > 1. This arises from the failure to consistently 
incorporate the elastic threshold in moments of structure 
functions at finite Q 2 , resulting in nonuniformity of the 



Q 



oo and n — > oo limits, where n is the rank of 



the moment. After performing an inverse Mellin trans- 
form on the moments, the extracted structure functions 
consequently acquire incorrect support at large xb [H|- 
A number of attempts have been made to redress the 
threshold problem by considering various prescriptions 
to tame the unphysical behavior as -> 
These approaches are not unique, however, and some- 
times introduce additional complications (see Ref. [l6[ 
for a review). 

An alternative approach, which avoids the threshold 
ambiguities from the outset, involves formulating TMCs 
directly in momentum spac e fl7l using the collinear fac- 
torization (CF) formalism jl8l. 19|. This method was 
heuristically applied by Aivazis, Olness & Tung [2(| and 
by Kretzer & Reno [21J to spin-averaged structure func- 
tions. More recently Accardi & Qiu [17| applied this for- 
malism to deep inelastic structure functions at large xb, 
carefully taking into account the elastic threshold and 
thereby solving the threshold problem. However, in the 
handbag approximation, without introducing a suitable 
jet function accounting for the invariant mass of the final 
hadronic state [13] , leading order structure functions can 
still be nonzero at xb = 1. 
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In this Letter we use the CF framework to derive target 
mass corrections to the leading twist g\ and g-i structure 
functions and the A\ polarization asymmetry. In Sec. Hi] 
we outline the main steps in the derivation; a more de- 
tailed account will be presented elsewhere [22:]. In Sec. IIIII 
we compare and contrast the predictions for the TMCs 
in CF with those using the standard OPE formulation. 
We find that the TMCs using the CF method are gener- 
ally larger for the g\ structure function than in the OPE. 
However, since the TMCs are qualitatively similar for gi 
and Ft, the effects largely cancel in the A\ asymmetry, 
although the residual effects can still be up to 20% at 
large xb, and for A\ even differ in sign for the CF and 
OPE approaches. Finally, in Sec. IIVI we summarize our 
findings and preview future work. 



II. TMC IN COLLINEAR FACTORIZATION 

The computation of TMCs in collinear factorization 
makes use of the factorization theorem relating the 
hadronic tensor for j*N scattering to the partonic 
tensor wj" for the scattering of a virtual photon from 
a parton of flavor /. The target mass corrected struc- 
ture functions are then obtained by suitable projections 
of the hadronic tensor without neglecting the target mass 
M relative to Q 2 at any stage. 

The hadronic tensor for spin-dependent inclusive scat- 
tering of leptons from nucleons is given by 



W^(p,q) = —s^q p 
p-q 



S a gi (xb,Q 2 



S ■ q 
1 

p ■ q 



Pa ) 92 (%B , Q ) 



(1) 



where p and q are the target nucleon and virtual photon 
four-momenta, respectively, and S is the nucleon spin 
vector, with S 2 = —M 2 and S ■ p = 0. We work in 
collinear frames, defined such that p and q do not have 
transverse momentum. This allows us to decompose p, q 
and the parton four-momentum k in terms of light-cone 
vectors and as HH 



fractional momentum 



£ = -Si = 2xb 

p + i + v/i + r* 



(3) 



coincides with the Nachtmann scaling variable Q, with 
7 2 = Ax 2 B M 2 /Q 2 . In the Bjorken limit (Q 2 ^>oo at fixed 
xb), £, — ► xb and we recover the standard kinematics in 
the M « approximation. 

The nucleon polarization vector S can be decomposed 
into longitudinal (Sl) and transverse (St) components, 



S" 



aX S£ 



Qfl 

>Jrp 



(4) 



where S\ = 



M , Sl ■ P — St • P — 0, and the nucleon 
hclicity A = ±1 indicates polarization parallel or antipar- 
allel to the nucleon direction of motion. The degree of 
longitudinal polarization is given by a = y/l + S T /M 2 



1 describing nucle- 
■M 2 ) or longitudinal 



with the limits er m i n = and cr max - 
ons with purely transverse (S T = - 
(S T — 0) polarization, respectively. 

Collinear factorization for the hadronic tensor can be 
obtained by expanding the parton momentum k around 



its on-shell (k 2 
poncnt, 
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0) and collinear (kr — > 0) com- 



k' 1 



k" 



xp n H 



(5) 



In terms of the on-shell parton momentum k we can then 
define the collinear invariant 



2k ■ q 



£ 

x 



(0) 



where the second equality hold for massless quarks, to 
which we restrict this analysis. 

The spin vector of a collinear parton, s M , is defined 
analogously such that s 2 — and s-k — 0. For a massless 
spin-1/2 quark, as well as for a massless spin-1 gluon, the 
spin vector can be written as 



A/ A W 



(7) 



p^ = p + n tJl + 



M 2 
2p+ 



2^ 



(2) 



fc" = xp+n» + k l +k T n » + K ■ 



2xp^ 



where n 2 



and n • 



1. The transverse parton 



momentum vector k!f satisfies kr ■ n = kT 



0. 



The nucleon plus- momentum p + = (po +p3)/v / 2 can 
be interpreted as a parameter for boosts along the z-axis, 
connecting the target rest frame to the hadron infinite- 
momentum frame. In terms of the plus-components of 
the momenta, the parton fractional light-cone momen- 
tum is defined as x = k + /p + , while the virtual photon 



where the parton helicity A/ = ±1 corresponds to a par- 
ton with spin parallel or antiparallel to the proton longi- 
tudinal spin. 

According to the QCD factorization theorem the 
hadronic tensor can be factorized as 



W^(p,q,S) 



E 



i/xB dx 



w 



f V {k,q,s)<p/ (x,Q 2 ) , 



(8) 

where w^ v is the partonic tensor for scattering from a 
parton of flavor /. (Note that the vector s in the ar- 
gument of w^ v depends on Xf.) The upper limit of in- 
tegration in Eq.©, viz., x max — £/xb, guarantees that 
structure functions vanish for xb > 1 [17] . in contrast 
to Refs. [13, [2l| where x max = 1. The neglected terms 
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of order 0(k^ — k^) in the collinear expansion are sup- 
pressed by powers of A 2 /Q 2 , with A some hadronic scale, 
and contribute to the restoration of gauge invariance in 
higher- twist diagrams [23| . The factorized expression © 
is obtained with the additional approximation of neglect- 
ing the intrinsic parton fcy and parton off-shcllness in the 
kinematics of the handbag diagram. The (x B , Q 2 ) region 
where this approximation is valid has been estimated in 
Ref. [ItJ ■ A detailed account of non-zero kx requires go- 
ing beyond the collinear factorization formalism used in 
the present analysis. 

Polarized scattering is described by the antisymmetric 
part of the tensor, which can be decomposed in terms of 
the partonic gij and g 2 j structure functions, 



wf(k,q) 



k ■ q 



e^ pa q P 



s <?9ij{xf,Q 2 ) + (s CT - ^-k a )g 2J (x f ,Q 2 ) 
v k ■ q ' 



(9) 



For ease of notation, in the following we will omit the 
dependence on Q 2 of the structure functions and parton 
distributions functions. The function ip^ f in Eq. ([5]) is 
the parton distribution function for a parton of flavor / 
and hclicity Xf inside a nucleon. In the light-cone gauge, 
and at leading order in a s , this is defined as 

(10) 

where ipf is the quark Dirac field. For polarized scatter- 
ing the spin-dependent quark distribution function Aipf 
is then given by 



1 



A<p f (x) = - [tpj (x)-<p f (x)] 



(11) 



Note that the factorized expression in Eq. ((Sj) is suit- 
able for discussing the contribution of helicity parton 
distributions to the hadronic tensor at leading order in 
the expansion of parton correlators in powers of l/p + 
[24l 12a , |26| . Extension to transversity distributions, 
or inclusion of higher order corrections in l/p + (corre- 
sponding to "dynamical twist" > 3 in the language of 
Refs. [2J, require a generalization of Eq. ( El). Even 
at 0(l/p + ) the TMCs can become nontrivial [lq . and 
we will discuss these higher order corrections elsewhere. 

Using suitable projection operators, structure func- 
tions can be projected from the hadronic and partonic 
tensors in Eqs. (fTJ) and ([9]), and using Eq. §8$ one finds 



9i 



g2(x B ) = -9i(xb) 



(12) 
(13) 



for the target mass corrected structure functions at 0(1) 
in l/p + . At leading order in a s the partonic structure 
function gij is proportional to 5(x — £), in which case 
the target mass corrected nucleon g\ structure function 
is given by 



9x(x B ) 



1+7 



— 2 9?\0 , 



(14) 



where g^ is the structure function in the massless target 



limit, M 2 /Q 2 ^0. Note that Eq. $T4]| is strictly valid only 
at leading order. At higher orders the massless limit g\ 
structure function generalizes to 



9?\xb) 



E 



1 dx 



T*'(-W(v)» (15) 

X \ X J V X J 



with <72°' (xb) = ~9\ > \ x b)- Clearly, in general one has 
9i{xb) 7^ 9i (0 because of the 1/(1 + 7 2 ) factor in 
Eq. p2[) . and the upper limits of integration (i.e., £/xb 
versus 1). 

In actual polarized deep-inelastic scattering experi- 
ments one typically measures not the structure functions 
directly, but the virtual photon polarization asymmetries 
A\ and A 2l defined as ratios of spin-dependent to spin- 
averaged structure functions, 



Ai{x B ) 
A 2 (x B ) 



gi(x B ) ~j 2 g 2 {x B ) 
Fi{x B ) 
gi(x B ) + g 2 (x B ) 
7 Fx{x B ) 



(16) 
(17) 



Using the results in Eqs. (fl2^) - ||13p one can write the 
asymmetries in collinear factorization as 



Ai{x B ) 
A 2 (x B ) = , 



(1+7 yi(* s ) 



(18) 
(19) 



where F\ is the spin-averaged structure function, which 
in collinear factorization is given by (l7| 



(20) 



using a shorthand notation for the integral over x as 
in Eq. (TT2]) . The function ip f is defined as the sum of the 
helicity distributions in Eq. (JTUJ) , ff(x) = cpi(x)+(pj(x). 

The polarization asymmetry in collinear factorization 
can then be written at 0(1) in the l/p + expansion as 



Ai(x B ) 



E/ (jhj g Ay/) (Q 

E/ (*!,/ ® v>f) (0 



(21) 



Note that the (1 + 7 2 ) prefactor is absent if the asymme- 
try is written in terms of the parton distributions directly. 
In the M 2 /Q 2 — > limit the asymmetries are given by 



AO), x _ gi (x B ) 
Fi '(xb) 



A 2 0) (x b )=0 



(22) 
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The massless asymmetry is also directly related 
to the lepton asymmetry Am for scattering leptons with 
longitudinal polarization aligned and anti aligned with 
the nucleon polarization, 



4 .> (IB) = 4!|£> 



(23) 



where D is a depolarization factor of the virtual photon 
[27| . A commonly used approximation in experimental 
data analysis relates the longitudinal lepton asymmetry 
with the ratio of the g\ and F± structure functions, 



expression, Eqs. ([T2"]) and (fT4| . The factor (I + 7 2 ) -1 can 
be traced back to the tensor decomposition of lf F and 
has the same origin as the factor appearing in Eq. (|12jl . 
while the remaining (1 + r y 2 )~ 1 l 2 XB / factor arises from 
the OPE treatment of TMCs. Substituting Eq. ([27]) in 
Eq. (fl6|) one obtains for the A\ asymmetry in the OPE: 



A? PB (x B ) 

(1+7 2 ) 



F° PE (x B ) 



9? PE (xb) 



y 



(28) 



(24) 



From Eqs. (fT8|) and ([23]) this is equivalent to assuming 
that 



A, » ^ 



(0) 



(25) 



In the next section we shall test the validity of this ap- 
proximation numerically, and compare the results of the 
collinear factorization with the target mass corrections 
obtained from the OPE. 



III. COMPARISON WITH THE OPE 

A common prescription for evaluating target mass cor- 
rections uses the operator product expansion to com- 
pute moments of structure functions at leading twist, 
including the trace terms which introduce the kinemati- 
cal M 2 /Q 2 corrections, and extracts the TMC structure 
functions through an inverse Mellin transform d, 0, H, 

mm- 

The resulting target mass corrected g\ and g 2 
structure functions can be written as [8j: 



^ PE (^) = 7jT^y^) 

7 2 [ 1 dvrx B +Z 7 2 -2 



(1 + 7 2 ) 2 A v 



t ' 2v /TT^ l0g( ? 



(26) 

u^-ru)- -nY^-ru) /' %? PE (y)- (27) 



x B y 



The expression for g$ m Eq. ([27]) is known as 
the Wandzura-Wilczek relation [28|, and was shown in 
Ref. Q to survive target mass corrections. This expres- 
sion differs from Eq. ([T3]) . obtained at 0(1) in the l/p + 
expansion in collinear factorization, by the presence of 
the integral term. In collinear factorization such a term 
emerges at 0(l/p + ), and Eq. (12"T|) holds if one neglects 
quark-gluon-quark correlators [29|, |30| and matrix ele- 
ments related to the Wilson line in the expansion of the 
quark-quark correlators (3lj |. 

The prefactors for in the first term of Eq. ([26]) dif- 
fer from those in the corresponding collinear factorization 



which again differs from Eq. (JT8J) in the integral term. 
One should also note that the Wandzura-Wilczek rela- 
tion (|27|) is not a direct consequence of the OPE (2f| H2] , 
which leaves open the possibility of 5-function contribu- 
tions at xb = to the right-hand-side of Eq. (f27|) . To 
explore the phenomenological consequences of contribu- 
tions to A\ from subleading powers of l/p + , we consider 
both definitions in the numerical evaluation of A\. 

In Fig.[T]we compare the results of a leading order eval- 
uation of target mass corrected versus uncorrected proton 
gi (left panel) and F\ (center panel) structure functions, 
and polarization asymmetries A\ (right panel) for several 
Q 2 values, using the leading order GRSV2000 (standard 
scenario) polarized parton distributions for g\ 1331 and 
the GRV98LO unpolarized distributions for F\ [34]]. For 
both the CF and OPE corrections, the g\ ratio dips below 
unity at intermediate xg, 0.2 < < 0.5, before rising 
dramatically at larger ig . The magnitude of the dip and 
the steepness of the rise for are naturally greater at lower 
Q 2 . However, while the size of the TMCs at xb % 0.5 is 
< 2-3% for Q 2 > 10 GeV 2 , at larger xb the corrections 
remain significant even at much larger Q 2 . For these 
reasons, the commonly adopted cut Q 2 > 1 GeV 2 for 
polarized parton distribution function analysis requires 
inclusion of TMCs for extracting precise PDFs. 

For the A\ polarization asymmetry the TMC effects 
largely cancel in the ratio because the TMCs in the 
F\ structure function are similar to those in (1 + 7 2 )gi 
S US [HI, Nevertheless, the residual effects can still 
be up to 20% at Q 2 = 1 GeV 2 , decreasing to - 2-3% 
at Q 2 = 10 GeV 2 . This provides a quantitative test 
of the validity of the commonly used approximation in 
Eq. (|24|) for the longitudinal asymmetry A\\ in terms of 
A\. Interestingly, the TMC effects drive A\ in opposite 
directions for < 0.7, with A\ increasing relative to 
Ai in the OPE approach but decreasing in the CF for- 
mulation. This is due to the fact that for most xb values 



> gi, while F^^ < F 1 CF . Such ordering arises 



mainly from the different prefactors for g\ in Eq. (|12[) and 
in the first term of Eq. [55] (for the analogous formulas for 
F x see Ref. [H). 

The effect of using the Wandzura-Wilczek relation in 
the computation of A\ instead of the CF result g\ = —52, 
indicated by the shaded band in Fig. [T] has < 5% effect 
in general, and is negligible for Q 2 > 3 GeV 2 . The contri- 
bution of the Wandzura-Wilczek term is small compared 
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FIG. 1: Ratio of target mass corrected to massless proton g\ (left panel) and F\ (center panel) structure functions and A\ 
polarization asymmetry (right panel) in collinear factorization (solid) and in the OPE (dashed), for Q 2 — 1, 3 and 10 GeV 2 . 
For Ai the shaded band for the OPE result indicates the effect of using the Wandzura-Wilczek relation, Eq. (|27[) (lower bound), 
or the identity gi + g2 = 0, Eq. (|13[) (upper bound). 



with the differences between the two TMC schemes con- 
sidered. 

The differences between the two TMC implementa- 
tions can be seen more dramatically in Fig. [2j where 
the ratios of OPE and CF target mass corrected g\ 
(left panel) and A\ (right panel) are presented. At 
Q 2 = 1 GeV 2 the gi structure function corrected using 
the OPE prescription can be up to ~ 20% larger than 
that using the CF approach at ~ 0.4, with the dif- 
ference decreasing at larger x b ■ The differences diminish 
with increasing Q 2 , so that by Q 2 — 10 GeV 2 the meth- 
ods give essentially the same results at the 2% level for 
all x B < 0.8. 

The polarization asymmetry is similarly found to be up 
to ~ 20 — 30% larger within the OPE approach at xb > 
0.7 for Q 2 — 1 GeV 2 , depending on the prescription used 
for g2, but again decreasing to < 2% for Q 2 = 10 GeV 2 . 
In all cases the TMCs are larger for Ai in the OPE than 



in the CF approach. These results clearly highlight the 
need for a careful treatment of TMCs in the low-Q 2 and 
large-a;^ kinematics. 



IV. CONCLUSION 

In this study we have derived for the first time the 
target mass corrections to the spin-dependent nucleon gi 
and gi structure functions, as well as to the polariza- 
tion asymmetry A±, in the framework of collinear fac- 
torization. In the CF framework the threshold problem 
affecting the OPE framework is naturally avoided by di- 
rectly implementing four-momentum conservation in the 
handbag diagram, rendering the structure functions zero 
for xb > 1. A further advantage of this formalism is 
that it can be readily extended to processes such as semi- 
inclusive DIS, where the OPE is not available, and indeed 
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FIG. 2: Ratio of the gi structure functions (left panel) and A\ polarization asymmetries (right panel) computed with target 
mass corrections in the OPE and CF formalisms at Q 2 = 1 (largest ratios), 5 and 10 GeV 2 (smallest ratios). For Ai in the 
OPE is as in Fig. Q] 



to any other hard scattering process. Additional correc- 
tions to structure functions at large xb, such as from jet 
mass corrections or threshold resummation, can also be 
naturally incorporated together with TMCs. 

The numerical results for the target mass corrections to 
the leading-order g\ structure function in CF are found to 
be qualitatively similar to those obtained from the OPE, 
but up to 20% larger at Q 2 = 1 GeV 2 . The corrections 
become smaller at larger Q 2 , with differences between the 
CF and OPE results < 2-3% for Q 2 = 10 GeV 2 . Nev- 
ertheless, the TMCs remain significant at xb > 0.7 even 
for Q 2 > 10 GeV 2 , and need to be taken into account 
when analyzing large- xb data. The numerical difference 
between the two schemes is likely to increase in a next-to- 
leading order computation, where the convolution over x 
in Eq. (j 1 2|) is performed only up to (,/xb instead of 1 
17]. 

Since the TMCs are qualitatively similar for the g\ 
and F\ structure functions, they largely cancel in the A\ 
asymmetry, although the sign of the correction is oppo- 
site in the CF and OPE approaches over most of the 
range of xb- The CF target mass effects in A\ can be 
as large as 20% at Q 2 = 1 GeV 2 for x B ~ 0.8-0.9, again 
decreasing to less than a few percent by Q 2 = 10 GeV 2 . 
The commonly used approximation relating A\ directly 
to the longitudinal asymmetry A\\ will therefore break 
down at low Q 2 , so that accurate determination of po- 
larized structure functions will require measurement of 
both Ay and the transverse asymmetry A± [27l[3^|. 

In the future, this analysis can be extended in several 
directions. Firstly, while the CF formalism avoids 
unphysical regions in dealing with the threshold prob- 



lem, the corrected structure functions remain nonzero 
at xb = 1. To tame this behavior one can follow the 
approach of Ref. by introducing jet mass corrections, 
which render the TMC structure functions zero in the 
limit ib -t 1. Furthermore, while we have restricted 
ourselves to massless quarks, the generalization to 
heavy flavors can be accommodated within the collinear 
factorization framework. In addition, future quantitative 
analysis of large-xs and low-Q 2 data will require TMCs 
to be computed for structure functions at subleading 
powers in l/p + , which will be necessary for a more 
complete treatment of 172, for instance. Finally, work 
is currently in progress (3o| to extend the collinear 
factorization formalism to semi-inclusive deep inelastic 
scattering, where the target and hadron mass effects are 
yet to be evaluated, and it will be interesting to address 
the case of transverse momentum dependent parton 
distributions (24[, which will enable the role of the 
parton intrinsic transverse momentum to be quantified. 
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